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1. Introduction 

Effective field theory (EFT) has proven to provide a useful tool for describing a wide 
class of meson-meson, meson-nucleon, and nucleon-nucleon (NN) processes with and 
without external probes in the low energy regime [1]. However, special care for the NN 
processes is required to deal with the 1 S'o channel for its long scattering length (or existence 
of a quasi bound state) and the 3 S'i — 3 Di channel for existence of the deuteron. In general, 
EFT is based on the perturbative expansion of physical observables in terms of small 
external momenta, but a non-perturbative treatment would be required for positions of 
the singularity at the small energy scales, which associated with the long scattering length 
and the small binding energy, compared to the chiral symmetry breaking scale, pion mass 
m w . 3 To deal with the problem, Weinberg has suggested counting rules that allow one 
to handle this non-perturbative problem and derive the NN potential systematically [3]. 
With the potential calculated to a given order, the S'-matrix is calculated from a wave 
function obtained by solving the Schrodinger equation. This scheme has shown good 
accuracy and convergence with only a few leading terms [4, 5, 6, 7]. Kaplan, Savage, and 
Wise (KSW) suggested new counting rules, in which pions are treated perturbatively and 
do not appear in leading order (LO). They also employed the power divergence subtraction 
scheme for regularization [8]. Only a LO contact two-nucleon interaction is treated non- 
perturbatively and thus physical observables can be directly calculated from Feynman 
diagrams expanded order by order. Over the last decade both the Weinberg and KSW 
schemes have been used extensively in studying few-nucleon systems. 4 For reviews, see, 
e.g., Refs. [10] and [11]. 

Convergence of deuteron observables becomes slow near the deuteron pole due to a 
large expansion parameter 7 pd ~ 0.4, where 7 ~ 45.7 MeV and pd ~ 1.764 fm in the 
KSW scheme. This would be a typical expansion parameter of the deuteron channel both 
in pionless effective field theory (pionless EFT) [12, 13], where the pions are regarded as 
heavy degrees of freedom in study of very low energy reactions and integrated out from 
effective lagrangian, and in theories that treat pions perturbatively. It was suggested, 
however, that adjusting the deuteron wave function to fit the asymptotic Estate nor- 
malization constant, Zd = 7 Pdj (1 — 7 Pd), assures more efficient convergence [14, 15, 16]. 
By introducing a dibaryon field which represents a resonance state (or a bound state) of 
two nucleons 5 , Beane and Savage showed that in terms of dibaryon EFT (dEFT) with- 
out pions, a long tail of the deuteron wave function can be naturally derived from the 
two-nucleon(dibaryon) propagator at the deuteron pole [19]. It was also shown that, with 
the use of dEFT, a good number of diagrams at a given order even in pionless EFT can 
be cast into a few diagrams. This feature is expected to be useful when higher order 

3 Even if the scattering lengths were of the order of pion mass, still nonperturbative treatment of the 
NN amplitude would be required. See Rcf. [2]. 

4 Detailed and careful studies of both of the schemes have been reported in Refs. [2, 9]. 

initially, the dibaryon field was introduced by Kaplan in the 1 So channel NN scattering [17]. Subse- 
quently, the dibaryon formalism without pions was applied to studies of three nucleon systems by Bedaquc 
and van Kolck[18]. It was shown that the pole structure of the amplitude at the low energies can be more 
efficiently reproduced with the use of the dibaryon field and resummation of effective range. 
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corrections need to be calculated. One can also expect that dEFT should be an attrac- 
tive approach in attempts to incorporate, in the framework of EFT, radiative corrections 
(photon degrees of freedom) for the two-nucleon reactions [20, 21]. 

In applying dEFT to reactions involving an external probe, there appear additional 
low energy constants (LEC's) in effective lagrangian. LEC's that enter into electromag- 
netic (EM) interaction on the two-nucleon system appear in the vector(photon)-dibaryon- 
dibaryon (Vdd) interactions in terms of the dibaryon fields and can be fixed, e.g., from 
the total cross section of radiative neutron capture on a proton, np — > dj, or from the 
deuteron EM multipole moments. 6 These LEC's are considered to subsume the high en- 
ergy physics, such as meson exchange currents or heavy meson exchanges, that has been 
integrated out from the effective lagrangian. Using Weinberg's counting rule with pionful 
theory, these contributions from the meson exchange currents are higher order corrections 
and are known that their corrections are numerically about a few ~ 10% in the reactions 
considered in this work. Similar results have been reported in the KSW scheme and pi- 
onless EFT (the references will be given in each section.) However, as we will show, the 
contributions from the LEC's of the Vdd vertices in dEFT turn out to be anomalously 
larger (about 40 ~ 50% of LO contributions) than the former EFT results. 

Each version of EFT has its own expansion scheme and values of LEC's for a given 
vertex can differ from one version to another. The large (40 ~ 50%) cotributions from 
the LEC terms in dEFT may be regarded as consistent with the expansion parameter 
(7p d ~ 0.4) in the pionless theories [12, 13]. However, such large contributions from the 
LEC terms with the EM probe could alter our knowledge about this class of the LEC's 
and imply that convergence can be much slower than that found in the previous EFT 
calculations. This could be a serious drawback in calculations with dEFT, since to obtain 
sufficient accuracy, one requires a higher order calculation due to the slow convergence. 
Subsequently, taking into account one higher order can be a formidable task and more 
importantly, the slow convergence of dEFT can make the uncertainty of the results at 
a given order larger than that in the other EFT calculations. In order to make the 
convergence of dEFT more efficient, we propose a simple way to handle this class of the 
LEC's, i.e. the LEC's of the Vdd vertices, in dEFT. 

In this work, we explore several applications of dEFT for the deuteron reactions. We 
consider a pionless effective lagrangian that includes the dibaryon fields in the 1 S'o and 
3 S\ — 3 Di channels. Those LEC's appearing in the strong interaction part of the dibaryon 
lagrangian are fixed by the effective range parameters in each channel. The EM interac- 
tions between the dibaryon fields and an external vector field introduce additional LEC's 
in both isoscalar and isovector transitions. At first, we fix the LEC's of the Vdd vertices 
by experimental data. The contributions from them, as mentioned before, turn out to 
be anomalously larger than those found in the previous EFT calculations. We suggest 
modified counting rules in which these LEC's are enhanced by an order of Q -1 . In the 
new counting scheme, we suppose that a prime portion of the LEC's of the Vdd vertices 
in dEFT originates mainly from one-body vector(photon)-nucleon-nucleon (VNN) inter- 
recently Detmold and Savage [22] suggested that lattice simulations allow us to estimate the LEC's 
that feature in dEFT describing reactions with an external probe. 
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actions. We find that this scheme for fixing the LEC's makes the convergence of dEFT 
as efficient as other EFT and well reproduces the deuteron EM form factors and the total 
cross section for np — > d'-f obtained in the effective range theory (ERT) [23]. Including 
higher order terms in our dEFT calculations, we find that the results of dEFT become 
comparable to those of potential model calculation using Argonne vl8 potential [24]. 

This paper is organized as follows: In Sect. 2 the effective lagrangian is given. We 
calculate in Sect. 3 the EM form factors of the deuteron and fix the LEC's in the EM in- 
teraction lagrangians. We then calculate physical observables in the elastic e-d scattering, 
and compare our results with those of ERT and a potential model calculation, and also 
with the experimental data. In Sect. 4, we investigate the np — > d'-f reaction. The final 
section, Sect. 5, is devoted to discussion and conclusions. In Appendix, we give details 
on the two-nucleon(dibaryon) propagators and NN scattering amplitudes, and determine 
the LEC's in the strong-interaction lagrangian. 

2. Effective lagrangian with dibaryon fields 

A pionless effective lagrangian for the nucleon and the dibaryon fields interacting with 
an external vector field can be written as 

C = C N + C s + C t + C st , (1) 

where £jv is the one-nucleon lagrangian, C s and £ t are the lagrangians for the dibaryon 
fields in the 1 S'o and 3 Si channels, respectively. C st is the lagrangian that accounts for 
the transition between 1 Sq and 3 Si channels through the isovector EM interaction. 7 
Cn in the heavy-baryon formalism reads 

(„ . D f — D 2 — SI ((1 + « v )/+ + (1 + K 3 )f%)_ 

+2d 7 v»[D», /£] + 2e 54 ^, S»][D\ [D x , /£]] + • • •} N , (2) 

where the ellipsis represents terms which we do not include in our calculation, v M is the 
velocity vector satisfying v 2 = 1; we take = (1, 0). is the spin operator 2S^ = (0, a). 

= <9 M — \t • V M — | Vf , where V M and are the external isovector and isoscalar vector 
currents, respectively. /+, = \ ■ {d^V u - d u V^ and f^ v = \ (<9 M Vf - d„Vjp) . m N is the 
nucleon mass and Ky (ks) is the isovector (isoscalar) anomalous magnetic moment of the 
nucleon; k v = 3.70589 (k s = —0.12019). The terms of d 7 and e 54 correspond to those of 
d 7 and in the pionful one-nucleon lagrangian [25] and the LEC's d 7 and are fixed 
by isoscalar radii of nucleon EM form factors. 

7 The counting rule we employ here is the same as that in Rcf. [19]: we do not expand the amplitudes in 
terms of the effective range, pd or tq. Introducing an expansion scale Q < A(~ m n ), we count magnitude 
of spatial part of the external and loop momenta, |p| and |Z|, as Q, and the time component of them, p° 
and I , as Q 2 . Thus the nucleon and dibaryon propagators are of Q~ 2 and a loop integral leads to Q 5 . 
The scattering lengths and effective ranges are counted as Q <~ {7, l/a , l/pd, l/ r o}- Orders of vertices 
and diagrams are easily obtained by counting the numbers of these factors. LEC's are counted as Q° 
initially. However we will argue that some of them (the Vdd vertices) can be enhanced by Q~ x . 
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C s , C t , and C st for the dibaryon fields [19] may read 

1 
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(5) 



iCVp, where 



where the covariant derivative for the dibaryon field is given by = 
V^ xt is the external vector field. 8 C is the charge operator of the dibaryon fields and 

C = 0, 1, 2 for the nn, np, pp channel, respectively (where we have put e = 1). B is the 
magnetic field given by B = V x \> ext . The sign factors, a s and a t , turn out to be — 1 
(see Appendix for details). A t (A s ) is the difference between the dibaryon mass m t (m s ) 
in the 3 Si ( 1 S'o) channel and the two-nucleon mass; m ttS = 2mjv + A t)S . p d and r are 
the effective ranges for the deuteron and 1 S scattering state, respectively. P- S ^ is the 
projection operator for the S = 3 Si or 1 S channel; 



P [ 

-*- i 



Tr(PPP^) = (6) 



where U{ (r ) is the spin (isospin) operator. The operators for the D-state read 



ij,xy 



@if = — i DjD, - D 



O 



(2) 



$ix &jy g &ij $xy i 

~\ (b X Dy Pf Sl) + Pf Sl) DxDy ~ D X Pf Sl) Dy ~ Dy D t 



(7) 
(8) 



The LEC's, y s and yt, represent the dibaryon- NN (dNN) couplings in the spin singlet 
and triple states, respectively, and they contribute to the two-nucleon loop diagram for the 
two-nucleon(dibaryon) propagator. They, as well as A S)t and a Stt , have been determined 

8 Since the external vector fields coupled with charges of proton and neutron are given by + ViJ) 

and j(V^ — V^), respectively, \>^ x p f or n p channel of the dibaryon field, for example, is obtained by a sum 
of the charge operators of neutron and proton, V^ xt (np) = \{V^ — V^) + ^(V^ + V^) = . Similarly, 
onehasV^*(pp)=V^ + ^andV^ t (nn) = V*-V3. 
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from the effective ranges in the 1 Sq scattering and the deuteron states. accounts 
for 3 Dx state mixture to the 3 Si state and is fixed by the asymptotic D-S ratio rj sd (~ 
0.0254). Details are given in Appendix. 

We consider LEC's with the external vector field in the literatures. Cq, Li, and L 2 are 
the LEC's for the Vdd vertices, and the coefficients, f / y/m n Pd and f/y / m/\Fro, of these 
terms in Eqs. (4) and (5) are introduced by a rule in Eq. (ff) in Ref. [19] to convert 
the two-nucleon field (NN) in pionless EFT to a dibaryon field (t or s) in dEFT. In this 
work we do not derive the most general effective lagrangian with dibaryons and without 
pions. Cq is determined from the electric quadrupole moment of the deuteron, and L\ 
and L2 are determined from the total cross section of the radiative neutron capture by 
the proton and the magnetic moment of the deuteron, respectively. We introduce a new 
LEC Cm, which can be fixed from the radius of the magnetic form factor of the deuteron. 
In general, one also has LEC's L[ and L' 2 for the vector(photon)-dibaryon- iViV (VdNN) 
vertices by using the rule mentioned above. Since they are higher dimensional terms, they 
could be considered to account for a higher order correction to the leading contributions 
of EM transitions. In the next two sections, we determine the EM LEC's listed above 
and proceed to calculate the EM observables in the elastic e-d scattering. 

3. EM form factors of the deuteron 

The EM form factors of the deuteron have been intensively studied within the Wein- 
berg's approach [15, 26], KSW scheme [27], and pionless EFT [28]. The electric form factor 
of the deuteron has also been calculated in dEFT without pions [19]. In this section we 
consider the deuteron EM form factors in dEFT without pions. 

A deuteron state \p, i) specified by momentum p and spin % satisfies the normalization 
condition (p',i\p,j) = (2ti) 3 5^ (p' — p)5ij. The nonrelativistic expansion of the matrix 
element of the electromagnetic current up to 0(Q 3 ) is given as [28] 



{p',i\Jem\pJ) 

(P,i\Jem\PJ) 



'E' + E" 



2m 



M<l)W + P) k + F M (qM qj - Sfa) 



(9) 



where p' = p + q and q — \q\. e is the electric charge, rrid the deuteron mass, and E (£") 
the energy of the deuteron in the initial (final) state, n is the space-time dimensions, 
4. Fc(q), Fm{q), and Fq(q) represent electric charge, magnetic dipole, and electric 



n 



quadrupole form factors, respectively. These dimensionless form factors defined in Eq. 
(9) are conventionally normalized as 



F c (0) = 1, ^F M (0) = fi M , -\f q (0) = n Q , 



(10) 



where ^m{= Pd 2^7) i s the magnetic moment of the deuteron and /iq is its electric 
quadrupole moment p d = 0.8574382284(94) [29] and p Q = 0.2859(3) fm 2 [30]. The 
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(a) 



(b) 



Figure 1: Diagrams for deuteron electric form factor Fc(q) at LO, O(Q ). The single line 
stands for the nucleon field, the double line for the dibaryon, and the wavy line for the 
external vector field. The vertex with a dot coupled with the external vector field in (a) is 
proportional to v^, while that with a dot in (b) is proportional to v^a t . The dNN vertex 
with a dot in (a) is proportional to y t . 



charge radius of the deuteron y (r^ h ) is defined by Fc{q) = 1 — \(r'l h )q 2 + 
empirical value is (r 2 ch ) = 2.1303(10) fm [31]. 



and its 



3.1 Electric form factor 

Two LO, O(Q ), diagrams in Fig. 1 contribute to the electric form factor of the 
deuteron, Fc(q). Three-point vertex functions calculated from Fig. 1 (a) and (b) read 



iT 



-iSijV^at arctan 



Pdq 



q_ 



(11) 
(12) 



where i and j denote the final and the initial spin state of the deuteron, respectively; p, is 
the Lorentz index for the current. Multiplying the three-point vertices with the normal- 
ization factor Z d of the deuteron wave function (Z d is derived in Eq. (51) in Appendix), 
one obtains the charge form factor of the deuteron at LO as 



and thus 



v»F%° 



■'j 



IPd 



1-lPd 



L (ij) +L (ij) 



4 f q\ 

arctan — — 1 

Pdq V 4 7/ 



(13) 



(14) 



This result is the same as that calculated with ERT [27]. We obtain the same charge 
radius (r c \) as ERT, {r 2 ch ) ER = = (1.985) 2 fm 2 . 

Higher order corrections to Fc(q) have been studied, e.g., in pionless EFT [28], but 
a slightly different counting rule is employed in it where one expands the form factors 
in terms of the effective range, i.e., counts 7 ~ Q, 1/pd ~ Q°, and 7p rf ~ Q. (See also 
footnote 7.) After expanding our result in terms of (jpd), one can easily compare it with 
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(a) 



(b) 



Figure 2: Diagrams contributing to deuteron magnetic form factor -Fm(?)- The leading 
order diagram (a) is of 0(Q), while the next-to-leading order diagram (b) is of 0(Q 2 ). 
The vertex with "X" in (a) is proportional to (1 +/cs)/(2mjv) while that with "X" in (b) 
is proportional to L 2 . See the caption for Fig. 1 as well. 



Eqs. (3.21,22,23) in Ref. [28] and find agreement with each other up to the subleading 
order. A difference in the higher order comes out of the isoscalar vector radius of nucleon 
current and relativistic correction. 

One can verify that F^°(q) in Eq. (14) satisfies the normalization condition given by 
Eq. (10). It would be worth noting that the expression of F^°(q) has no free parameter 
because the leading Vdd vertex that contributes to the electric form factor stems from 
the covariant derivative of the dibaryon field (Fig. 1 (b)), which is proportional to the 
dibaryon charge C = 1 for the np channel and the overall factor o t . In other words, the 
magnitude of strength of the Vdd vertex is determined solely by the charges of nucleons 
associated with the one-body VNN vertex. At q — 0, the contribution of the diagram 
with Vdd vertex (Fig. 1 (b)) amounts to —40% of that with the VNN vertex (Fig. 1 
(a)). As will be shown in the next subsections, we find similar values of the Vdd to VNN 
ratios for the other form factors. More importantly, the quantity in the bracket of Eq. (14) 
exactly cancels the factor jp d /(l — jpa), the deuteron normalization factor Z d , at q = 0, 
which leads us to assign a new role to the LEC's of the Vdd vertices. 

3.2 Magnetic form factor 

Diagrams contributing to the magnetic form factor of the deuteron are depicted in Fig. 
2. The order of Fig. 2 (a) is LO (0(Q)) and that of Fig. 2 (b) is next-to-leading order 
(NLO) (0(Q 2 )) for the vertex functions. We will discuss however that the contribution of 
Fig. 2 (b) can be separated into LO and sub-leading parts, and show that this re-ordering 
makes the convergence more efficient. The three-point vertex for each diagram is 

iT§ )k = i(5 ikqj -5 jkqi )^^^cUn(j-y (15) 
where k is the index for the spatial part of the vector current. Multiplying the three-point 
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functions with the normalization factor Z d , one has the magnetic form factor Fm{q), 



1 TP f \ 7Pd 
F M (q) = 

2m d 1 - 7p d 



l + 4 / q \ 2L 2 

arctan — H 

2m N p d q XA'j J m N p d 



(17) 



At q = 0, 



2m, 



-i^(0) 



7Pd 



1 ~lPd 



1 + K S 1 2L 2 



(18) 



If one fixes the LEC L 2 using the experimental value of the deuteron magnetic moment 
Hm in Eq. (10), one obtains L 2 = —0.4033 fm. Comparing the magnitude of Fig. 2 
(a) with that of Fig. 2 (b), one finds that the latter is about —42% of the former. As 
mentioned in the previous subsection, this ratio is similar to the Vdd to VNN ratio in 
Fq°(0). Since the former EFT calculations [28, 32] show that the LO one-body term is 
dominant in the estimations of pm and the higher order terms give only a few percents 
correction, this rather large ratio we obtained here causes a difficulty to interpret L 2 in 
dEFT as the high energy contributions integrated out from the effective lagrangian. So 
we may divide L 2 into the sum of the LO and sub-leading contributions. 9 This modified 
counting of the LEC L 2 reads 



L 2 — Z/n -|- 5L 2 



L 2 



-^(1 + K S )p d , 



(19) 



where L\ is LO, while 5L 2 accounts for the sub-leading contributions. Inserting V\ into 
L 2 of the lagrangian, one finds that the coefficient of the Vdd interaction becomes (1 + 
K S )/{2m N ) which is the same as that of the one-body interaction. Note that L\ is of Q~ x 
because of the factor p^, so the diagram Fig. 2 (b) with L\ becomes the same order as 
the diagram Fig. 2 (a). With the modified counting of L 2 the magnetic form factor at LO 
becomes 



2m d 



F M (q) 



1 + k s -yp d 



L 2 =L° 



2m 



N 



1 + K S ^LO 
r, 



2m 



N 



IPd 
C (<?), 



— arctan ( — 
Pdq V 4 7, 



- 1 



(20) 



which is equivalent to the relation of -Fm(<z) an d Fc{q) obtained from ERT [27]. Numer- 
ically, one finds L\ ~ —0.3880 fm. This value is larger than L 2 fixed from Eq. (10) by 
about 3.7 %. This 3.7 % difference can be attributed to the contribution of 5L 2 (or the 
higher order term of L' 2 . We will discuss it below.) It is also easy to compare our LO 
result expanding in terms of ^p d with those of pionless EFT calculation, Eqs. (3.32,33) in 
Ref. [28]. We find good agreement with each other and that our SL 2 plays a similar role 
to the LEC in pionless EFT. Thus the partition of LEC L 2 in Eq. (19) is essential to find 
a relation between EFT's with and without dibaryon fields. 



3 The same partition has been employed for the LEC L 2 in the recent dEFT calculations [22, 33]. 
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(c) (d) 

Figure 3: Diagrams contributing to the deuteron magnetic form factor Fm{q) arising from 
the L 2 term. The "X" in the diagrams denotes the VdNN vertex corresponding to L 2 . 
Both diagrams belong to 0(Q 3 ). 



Next we consider some of higher order contributions from L 2 of the VdNN vertex 
(Fig. 3) and the D-wave of the deuteron (Fig. 4). Since there are no corrections to the 
isoscalar magnetic VNN vertex from the third order heavy-baryon chiral lagrangian [25] , 
we do not have NLO corrections to F M (q) from the one-body sector. 10 From the diagrams 
in Fig. 3 one obtains a three-point vertex function 



i(c+d) 
(ij)k 



iT^t,' = i(5 ik qj 



Sjkqi) 



27 



m N p d 



2ir 



2- 



(21) 



where dimensional regularization has been used for the loop integration. We find that 
contributions from the LEC's L' 2 and 5L 2 have the same momentum dependence in F M (q) 
and thus the two LEC's are redundant constants. In this work, we choose 5L 2 = and 
L 2 is fixed by the deuteron magnetic moment below. Sum of the diagrams in Fig. 4 gives 



?T (e+/) 
11 (<i)fc 



i(Sikqj - S jk qi) 



iT 



(a) 

(ij)k 



y/2(l + n s )r] sd 
32m N p d j 2 q 

i(Sikqj ~ S jk qi) (p' 2 + f) ■ 

m N p d 



A^yq — (I67 2 + 3g 2 )arctan 



47 



(22) 



Summing up the contributions from Fig. 2 to Fig. 4, we obtain the magnetic form 
factor of the deuteron with the NNLO corrections considered here as 



1 



2m d 



F M (q) 



IPd 



1-lPd 

16j 2 p d q 
2 7 
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2m jv 



-arctan 



Pdq 
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m N p d 



2tt 



L' 



(23) 



10 At NNLO, there are additional corrections from the radius and relativistic corrections to the isoscalar 
magnetic nucleon current. Since these corrections from one-body nucleon form factors are well known 
in the calculations of heavy-baryon chiral perturbation theory (see, e.g., Ref [34]), we will include the 
nucleon radii in a later subsection where the elastic e-d scattering is considered. 
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(e) (f) (g) 

Figure 4: Diagrams contributing to the deuteron magnetic form factor Fm(q) arising from 
the deuteron D-w&ve. The diagrams (e) and (f) are of 0(Q 3 ) and the diagram (g) is of 
0(Q 4 ), which will be considered later in e-d scattering. The "X" denotes the isoscalar 
magnetic vertex of nucleon-photon coupling proportional to (l+«s)/(2mjv), and the filled 
boxes stand for the deuteron D-w&ve in the diagrams (e) and (f). The filled diamond in 
the diagram (g) is a vertex proportional to Cm- 




(a) (b) (c) 



Figure 5: Diagrams for deuteron electric quadrupole form factor Fq(q). Order of the 
diagrams (a,b) is of 0(Q 2 ) and that of the diagram (c) is of 0(Q 3 ). The filled boxes in 
the diagrams (a,b) denote the D-w&ve contribution and the empty box in the diagram (c) 
represents a vertex proportional to Cq. 



where we have put 5L 2 = (i.e., L 2 = as mentioned before. Fixing L' 2 from Eq. (10) 
at q — 0, we obtain L' 2 = -0.07096 fm 5 / 2 , which gives about 2.5% correction to the leading 
term. We will study below a role of the higher order corrections to the form factor F M (q) 
in elastic e-d scattering. 

3.3 Electric quadrupole form factor 

Diagrams for the electric quadrupole form factor of the deuteron are depicted in Fig. 5. 
The diagrams (a) and (b) in Fig. 5 are of 0(Q 2 ), while the diagram (c) in Fig. 5 is one 
order higher, of 0(Q 3 ). As before, we will see the similar VNN and Vdd ratio in Fq(q) 
to that found in Fu(q)- Thus we employ a similar assumption to that is used for L 2 in 
Eq. (19) to extract the LO part of LEC Cq in the Vdd vertex in Fig. 5 (c). 
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Three-point vertex functions from the diagrams in Fig. 5 are obtained as 



Vsd 



(c)M 
(ij) 



4 v / 2Pd7V 



-A'jq + (I67 2 + 3g 2 )arctan 



m N p d 

Multiplying the three-point functions with the factor Z d , we obtain 
F Q {q) 




2m 2 



Pdl f_3 Vsd 

l-Pdl Uv^ 



pdl 2 q 3 



-A'yq + (I67 + 3q )arctan 



47 



C Q ) 



f k ,( 26 ) 

m N p d J 



and at q = 



Fq(0) _ 2p dl 



m„ 



I- Pdl 



1 Vsd Cq 



2 v / 2Pd7 3 m N p d 



(27) 



The normalization condition Eq. (10) leads to Cq = —1.709 fm 2 . Comparing the first 
and second terms in the bracket in Eq. (27), one finds that the latter is about —50% of 
the former. In addition, it shows a pattern similar to what we have seen in the cases of 
Fc(q) and FM{q)- One may notice, however, that the contribution to Fg(q) is generated 
by the D wave in the diagrams (a) and (b) in Fig. 5 and there is no one-body VNN vertex 
corresponding to the Vdd vertex of the LEC Cq. We assume here that the prime portion 
of the LEC Cq in dEFT, denoted by Cq, and the contributions from the two-nucleon 
loop diagrams (Fig. 5 (a) and (b)) are combined and produce a factor (1 — pdi) so that 
it cancels the factor (1 — Pdl)' 1 in Z d - We choose a modified counting of Cq as 



C 



Q 



c° Q 



5C 



Q- 



2s/2 7 



Vsd 

jm N p d , 



(28) 



where numerically C Q = —1.408 fm 2 . The contribution of 5Cq to pq, which can be 
attributed to the contribution of the high energy physics, is about —17.2%. This rear- 
rangement of Cq leads to more efficient convergence in dEFT calculation. 

Expanding Eq. (26) in terms of {iPd), we compare it with a result of pionless EFT, 
Eq. (29) in Ref. [35]. 11 Up to next-to leading order in the (iPd) expansion, we find good 

11 Two expressions of Fq(q) have been reported by the same authors in former pionless EFT calculations, 
Refs. [28] and [35], but these two calculations showed a slightly different result with each other. In Ref. [28] 
a LEC corresponding to C^"^ was fixed by using a mixing parameter ei in the 3 Si- 3 Di channel NN 
scattering, while in Ref. [35] it was fixed by using the asymptotic D-S ratio rj s d, which is the same quantity 
as what we have used for fixing the LEC ■ At q = 0, consequently, numerical results of leading order 
contributions to the deuteron quadrupole moment, fig , are different and obtained as /Jq° = 0.273 and 
0.335 fm 2 in Refs. [28] and [35], respectively. Moreover, a subleading order correction from /Jq° was 
obtained as ^(7Pd)A*Q° m Eq. (3.42) in Ref. [28]. The factor 1/2 of it cannot be obtained by expanding a 
factor 1/(1 — jpd) in Z d in terms of ("fPd) and appeared as another difference from the result in Ref. [35]. 
Therefore, we compare our result with that in Ref. [35] . 
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agreement with each other. In a higher order, a difference comes from nucleon radius 
terms which are higher order ones in our counting rule. This agreement here, however, is 
found without introducing the partition of the LEC Cq in Eq. (28). Since, as mentioned 
before, there is no one-body VNN vertex corresponds to the Vdd vertex of Cq, it may 
not be essential to introduce the partition in Eq. (28) for the LEC Cq. 

3.4 Elastic electron-deuteron scattering 

With the deuteron form factors obtained above we now make a brief study of elastic 
electron-deuteron scattering, e + d — > e + d. The differential cross section of the reaction 
is given by 



da 

dn 



da 

dn 



Mott I 



A(q) + B(q)tan 2 



(29) 



The form factors A(q) and B(q) are related to the deuteron form factors via 



A(q) = F 2 (q) + - V F 2 M (q) + - V 2 F 2 



Q 



(q), B(q) = ^ V (l+ V )F 2 M (q), 



(30) 



with i] = q 2 1 (4m 2 ,). 

In incorporating effect of the nucleon radius into the electric form factor F c (q), which 
have not been considered explicitly in the previous subsections, we introduce the radius of 
the isoscalar vector nucleon current and a radius of the dibaryon field into the diagrams 
(a) and (b) in Fig. 1, respectively. As discussed so far, order of a LEC of Vdd vertex, here 
for the dibaryon radius, should be modified and parted into LO and subleading part and 
the large potion of the LEC can be fixed by the nucleon radius. From the large portion 
of the LEC and the nucleon radius, in consequence, one has an expression as introducing 
the nucleon form factor (where we retain only the radius term) into F^°(q) in Eq. (14). 
From the small portion of the LEC one has a parameter to fix. Thus one has 



F c (q) 



Pdl 



1 - Pdl 



-± arctan ± 

PdQ \4~l i 



- 1 



(i-^(4> 2 )-^ 2 }, (31) 



where (r 2 Es ) is the isoscalar-electric radius of the nucleon Sachs form factor, 12 (r% s ) = 
0.604 fm 2 , and relates to the LEC dj in the one-nucleon lagrangian C N in Eq. (2) via 
( r Es) = — 6c?7- 5 C is a subleading radius correction from the dibaryon field. The corrected 



charge radius reads (r 2 h ) 



+ 



) + Za8 c . Using the experimental value of (r 2 h ), 



we obtain 5 C = 
to 



chl i~ Pdl 8 7 2 

0.011 fm 2 . Similarly, introduction of the nucleon radius in F M (q) leads 



2m, 



■F M (q) 



erf Pi 



1 - IPd I 2m 



N 



1 - 



6 



Ms 



PdQ 



-arctan 




- 1 



12 We have obtained the values of the nucleon radii (rj 



radii for the Dirac and Pauli form factors, in Table 1 of Ref. [36]. 



and (r 2 Ms ) from r[ 8 ^ and the isoscalar 
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^o Sd f 4 7g - (I67 2 + 3g 2 )arctan ( 



27 fmW Tl t 25C M 2 



+^J^L> 2 + ^ q * , (32) 

m N p d V 27T mArprf J 

where we have separated the LEC Cm = C m +5Cm (and one has C M = (l+Ks)pd('r Ms ) /24) 
and fixed the value of by the one-body radius, isoscalar-magnetic radius of the Sachs 
form factor of the nucleon, (r Ms ) = 0.598 fm 2 . In Eq. (32) we also have multiplied the 
.D-wave corrections from the diagrams (e) and (f) in Fig. 4 with the nucleon form factor 
because the photon couples directly with the nucleon in these diagrams. Since both of 
the corrections from the D waves and nucleon radius are of NNLO, this multiplication 
affects to terms in N 4 LO. (r Ms ) relates to the LEC e 54 in the one- nucleon lagrangian C N 
in Eq. (2) via (r 2 Ms ) = 24m N e u /(l + k s ). 

LEC 5Cm can be fixed from the magnetic radius of the deuteron, \J (r Md ) of F M (q). 
From Eq. (32) one obtains the relation for (r Md ), 



(r 



Md) 




Vsiisd , 2 



(1 - 7p d )87 2 



+ (4) 



12 7 <5C 



M 



[I - ^p d )m N 



(33) 



Precise empirical value of (r 2 Md ) is not available yet, so we choose a value of 5Cm, equiv- 
alently, a value of J (r Md ), as \J (r Md ) = 2.135 fm, so as to get a better fitting for experi- 
mental data of B(q) [38] up to momentum transfer 400 MeV. We ignore the radius effects 
in Fq(q) because the contribution of Fq(q) to A(q) is negligible due to the factor rf (see 
Eq. (30)) in the momentum range under consideration. 

The results for A(q) and B(q) are plotted up to q = 400 MeV in Figs. 6 and 7, 
respectively. Since our formalism is a pionless theory, the momentum transfers up to 
q = 400 MeV may be too large (compared to the pion mass) and it could not be suitable 
for the theory without pions. However, we show the higher momentum transfer region 
to compare our result with the existing low-energy experimental data of B(q). In both 
figures, the solid lines represent our present results, the dashed lines correspond to the 
results of ERT, and the dotted lines to those of a potential model calculation. In the 
potential model calculation, we employed the one-body (IA) current operator [37], the 
deuteron wave functions obtained from the Argonne vl8 potential (Avl8), and the nucleon 
radii from Ref. [36]. The experimental data are taken from Ref. [38]. Our LO dEFT 
calculations reproduce well the results of ERT in the small momentum region. With the 
higher order corrections included, the results of the dEFT approach become better than 
those of ERT and are comparable to those of the accurate potential model calculations 
for both A(q) and B(q) up to about q = 200 MeV. Values of r ch , p d and r M d calculated 
from ERT and Avl8 are compared with the ones used in our calculation in Table 1. 
Experimental values of them can be found e.g. in Table 1 of Ref. [39]. We should note 
that the experimental value of rud still has relatively large uncertainty compared to the 
accurate empirical value of r ch . More accurate measurements of B(q) are expected to 
reduce the uncertainty of rMd- 
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Figure 6: Form factor A(q) in the elastic e-d scattering. The solid line corresponds 
to our result, the dashed line to the result of ERT, and the dotted line to that of the 
potential model calculation that employs the one-body (IA) operator [37], the deuteron 
wave functions obtained from the Argonne vl8 potential, and the nucleon radii from 
Ref. [36]. The experimental data are taken from Ref. [38]. 




15 



CQ 




150 200 250 300 350 400 
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Figure 7: Form factor B(q) in the elastic e-d scattering. For further explanations, see the 
caption of Fig. 6. 



4. Radiative neutron capture by the proton 

Radiative neutron capture by the proton, n + p — > d + 7, has been intensively stud- 
ied within the frameworks of EFT employing the Weinberg's scheme [40] and the KSW 
scheme [41] and pionless EFT [28]; see also [19]. We study this reaction in dEFT in more 
detail here. 

Diagrams for the np — >■ d^ reaction are depicted in Fig. 8. The diagrams (a) and (b) 
in Fig. 8 are of LO, 0(Q 1/2 ), and the diagram (c) in Fig. 8 is of 0(Q 3/2 ), thus of a higher 
order. However, we will discuss below that, because of the Vdd vertex involved, it would 
be appropriate to extract LO contribution from the (c) diagram, which is of the same 
order as that of the (a) and (b) diagrams. Summing up the contributions of the diagrams 
(a), (b) and (c) in Fig. 8, one obtains the amplitudes for the initial 3 Si and states as 



27T7 



1 



: (7) 



iA (a+b+c) CS ) = 



1 - 7p d m N -n/-i p + |p d (7 2 + p 2 ) 



X 



l + «g l , 2 + 2n { 7 2 +P 2 L ' 

27717V 2 m N 



(34) 




2tt7 



1 



1 - 7p d m N -ip+ \r a p 2 
J + 7; r oP + -7, L i 

a 2 / 2777at 



(35) 
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ERT 


IA(Avl8) 


dEFT 


Exp. [39] 


Teh 


1.985 


2.115 


2.130* 


2.130 


rud 


1.985 


2.090 


2.135* 






0.880 


0.847 


0.857* 


0.857 



Table 1: Comparison of the physical observables; radii, r ch and r Md , of the electric and 
magnetic form factors and the magnetic moment fid of the deuteron. r c h and rud are 
given in units of fm. In the column representing dEFT, the entries with V are input used 
for fixing the LEC's. Experimental values can be found, e.g., in Table 1 of Ref. [39]. 




(a) (b) (c) 

Figure 8: Diagrams contributing to np — > dy. Diagrams (a) and (b) are of 0(Q 1 ^ 2 ), while 
diagram (c) is of C(Q 3/2 ). The single lines, the double lines with a filed circle, and the 
wavy lines denote nucleons, dibaryons, and photons, respectively. The double line with a 
filed circle, which is obtained in Fig. 10 in Appendix, represents the propagation of 3 Si 
and 1 S'o states in the intermediate state and the wave function normalization factor \fZ~d 
in the final state. 



Here p is the relative momentum of the two-nucleon system, and k is the momentum of the 
out-going photon; p — \p\, k — \k\, and k — kjk. e*, d \ and n are the polarization vectors 
for the out-going deuteron and photon, respectively, and q is the spin S = 1 vector for the 
initial 3 Si state. a and r are, respectively, the scattering length and effective range of 
the neutron-proton 1 S state. The LEC L 2 has already been fixed in the previous section 
using the deuteron magnetic moment fiM- The LEC L\ should be fixed by experiment. 

It is worth noting that the LO (0(Q^ 2 )) amplitude for the 3 Si state in Eq. (34) 
vanishes with the modified counting of L 2 , 

A (a +b +c) { 3 Si) =0 , (36) 

L 2 —L 2 

where the expression of L\ is given in Eq. (19). It is well known that the leading isoscalar 
transition between the scattering and bound states vanishes due to the orthogonality of 
the wave functions. We confirm that our treatment of the LEC L\ satisfies this condition. 
Higher-order terms can give non-vanishing contribution to the isoscalar amplitude, but 
since they are very small compared to the dominant isovector one [42], we neglect them 
in following calculations. 
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We determine the value of Li in Eq. (35) from the cross section of np — > d'-f at the 
thermal energy, a exp = 334.2 ±0.5 mb [43]. The total cross section in the CM frame reads 

,= ^^Y.\A\ (37) 



4p 



spin 



where a is the fine structure constant, and the CM energy E corresponding to the thermal 
neutron experiment is E — p 2 /m^ = 1.264 x 10~ 8 MeV. Using the formula for the spin 
summation J2 sp in\^(d) ' x ^(i))! 2 = ^ an d including only the amplitude for the 1 S 
channel, -4(%) in Eq. (35), we obtain L x = -4.427 ± 0.015 fm. 

It is to be noted that, if L\ is set equal to zero, i.e., L\ = 0, the relevant cross 
section would become a(Li = 0) = 502.3 mb, which is about 1.5 times larger than 
the experimental value, i.e., o{L\ = 0)/a exp ~ 1.50. The magnitude of the LO cross 
section in the KSW scheme (Eq. (36) in Ref. [41]) is smaller than o exp by about 13 %. 
Furthermore, the leading one-body operators in the potential model calculations also lead 
to a cross section that is smaller than a exp by about 10 %. As is well known, most of 
this 10 % deficiency can be accounted for by the meson exchange currents. Thus, here 
again we are facing the situation that, whereas the conventional treatments indicate the 
dominance of the LO contributions, the dEFT results exhibit uncomfortably large higher 
order corrections. To solve this problem, we assume that L\ is dominated by a leading 
contribution denoted by L\ which is chosen so as to reproduce the result of ERT, and that 
the rest, 5L\ contains information about the high energy physics that has been integrated 
out. Thus we consider the decomposition 

Li — L° + 5L U L\ = ~{\ + K V )(r + p d ), (38) 

where numerically L\ = —5.275 fm. If we insert the expression for E\ into L\ in the 
lagrangian, the coefficient of the Vdd term becomes - 1+Ky i»±fe which is the coefficient 
of the one-body VNN interaction apart from a factor (r + p d ) / (2y/r p d ) ~ 1.024 (which 
is very close to unity). Since L\ is 0{Q^ 1 ) due to the factor (r + p d ) in Eq. (38), the 
order of diagram (c) in Fig. 8 becomes the same as that of diagrams (a) and (b). At LO, 
the total cross section reads 



na(l + Kyf^al 

a{ l} ~ pm%{l- 1Pd ) 



1 - — - j(r + pd) 



(39) 



which is the same expression as that of ERT. Numerically we have o"^?) = 304.9 mb. 
This value is very close to a IA = 304.5 mb obtained in the potential model with the use 
of the one-body operator. 

Now we study higher order corrections for the np — > d'-f reaction. Higher order diagrams 
(0(Q 5 / 2 )) with the VdNN vertex proportional to L[ are depicted in Fig. 9. One may 
notice that the diagrams with the D-wave contribution vanish for the initial S-wave states 
and the corrections from the radii are negligible because of tiny momentum transfer in 
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Figure 9: Higher order diagrams for the reaction np — > d*y involving the VdNN vertex 
proportional to L! x . The order of diagrams is of C(Q 5/2 ). See the caption of Fig. 8 as well. 



the reaction. The combined contributions of diagrams (d,e,f) in Fig. 9 are given by 



Comparing the amplitude in Eq. (40) with that in Eq. (35), one finds that the momentum 
dependences of the amplitudes for the LEC's L\ and L\ are slightly different. This 
difference, however, does not affect the numerical results at small energies significantly. 
Here we compute the Ml transition contribution to the np — > d'-f cross section, using two 
sets of the parameters: 5L\ = 0.736 fm and L[ = 0, and 8L\ = and L[ = 4.464 fm 5//2 ; 
these LEC's have been fixed from a exp . Predictions for the np — > d'-f cross section at 
energies up to 1 MeV above the thermal neutron energy are given in Table 2. A fairly 
good agreement is seen between the results corresponding to the two different choices of 
the parameter set. Our dEFT results also show reasonable agreement with the recent 
EFT calculations [44, 45], as well as with the potential model calculation including the 
exchange current [46]. 

5. Discussion and conclusions 

In this paper we have reexamined dEFT without pions for the deuteron reactions 
involving the EM probe. The EM form factors of the deuteron and the total cross sections 
of the np — > g?7 reaction have been calculated. We have found that the LEC's of the Vdd 
vertices in dEFT give much larger contributions than the corresponding terms in the 
previous EFT approaches. Although it has not been completely explored why the LEC 
terms of the Vdd vertex give such large contributions in dEFT, we have proposed a 
practical prescription to extract the LO part of the LEC's in such manner that this LO 
part mimics the leading one-body VNN vertices. After the LO part is disentangled, the 
remainder is assumed to represent the high energy physics that has been integrated out. 
If one is to calculate it explicitly, one should take into account relativistic corrections, 
meson-exchange currents, higher partial waves, etc. However, an easy way to fix the 
higher order part of the LEC's is to fit it to experimental data. We were able to confirm 




(40) 
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E(MeV) 


(Ta/M \8La ) 
^ lvi i v i / 


(T A//"1 ( I 


a M i(EFT)\45] 


(TMifAvl8)[46l 

^ _M L \ a l / L J 


1.264 x 1(T 8 


334.2* 


334.2* 


334.2* 


334.2* 


5 x 10" 4 


1.668 


1.668 


1.667 


1.668 


1 x 10~ 3 


1.171 


1.171 


1.170 


1.171 


5 x 10" 3 


0.4953 


0.4954 


0.4950 


0.4954 


1 x 1(T 2 


0.3281 


0.3281 


0.3279 


0.3281 


5 x 10~ 2 


0.09814 


0.09820 


0.09810 


0.09820 


O.lOO 


0.04969 


0.04976 


0.04973 


0.04975 


0.500 


0.00775 


0.00782 


0.00787 


0.00781 


1.00 


0.0035 


0.0036 


0.0036 


0.00355 



Table 2: Ml transition contribution to the cross section for the np — > d'y reaction. The 
second and third columns are our dEFT results using the two sets of LEC's described 
in the text. The fourth column shows the results of the NNLO EFT calculation [45], 
while the fifth column gives the results of an accurate potential model calculation (Avl8) 
including the exchange current [46]. The LEC's in both dEFT and EFT and the strength 
of the exchange current in the potential model calculation are fixed at the values marked 
by *. 



that the higher order corrections defined in this way are indeed small, in conformity with 
the general tenet of EFT. 

With the LEC's thus determined, we have calculated the form factors A(q) and B(q) 
for elastic e-d scattering and found that the results of A(q) are very close to the exper- 
imental data for momenta up to 200 MeV. Furthermore, our results of A(q) and B(q) 
are comparable to those of the accurate potential model calculation. Our estimations of 
the total cross sections of the np — > d^ reaction for energies up to 1 MeV also agree well 
with the results of the other EFT calculations and the accurate potential model. With 
the proper treatment of the LEC's in the Vdd vertices, the convergence of dEFT becomes 
similarly efficient as other EFT's and this can be interpreted as an indication that dEFT 
is a useful tool for understanding a wide class of low-energy phenomena involving the 
deuteron and an external probe. 
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Figure 10: Diagrams for the two-nucleon(dibaryon) propagator at leading order: a double 
(single) line stands for the dibaryon (nucleon) field. 




(a) (b) 

Figure 11: Diagrams for the NN scattering amplitudes. Diagram (a) is for S- to S-w&ve 
channel and diagram (b) is for D- to S-wave channel. The dNN vertex with (without) 
a filled box denotes a vertex proportional to (yt,s)- The propagator of the two- 

nucleon(dibaryon) filed (a double line with a filled circle) is obtained from the diagrams 
in Fig. 10. 



In this appendix we rederive the two-nucleon(dibaryon) propagators and fix the LEC's 
using the known properties of NN systems in the S-wave channels; see also Ref. [19]. 

LO diagrams for the two-nucleon(dibaryon) propagators in the S'-wave channels are 
depicted in Fig. 10. Since the insertion of the two-nucleon one-loop diagram does not alter 
the order of the diagram [8, 12, 13, 18], the two-nucleon bubbles in the propagators should 
be summed up to infinite order. Thus the inverse two-nucleon (dibaryon) propagators in 
the center-of-mass (CM) frame read 



= i 



rriNVt, 
An 



2 r 



m N yL 



4:71 



Ana t<s E 

m N yl s 



ip 



(41) 



where \i is the renormalization scale of the PDS scheme, p is the magnitude of the nucleon 
momentum in the CM frame, and E is the total energy E ~ p 2 jm^. 

The two-nucleon (dibaryon) propagator in the spin singlet state, D s (p), can be renor- 
malized using the 1 S'o channel NN scattering amplitude. 13 The amplitude obtained from 
Fig. 11 (a) reads 



iA s = (-iy s )(iD s (p))(-iy s ) 



m N 



4tT(T s A s 

m N y'j 



ip 



(42) 



13 Although it is known that the expansion series of the ERT parameters in the 1 Sq channel converges 
well, we employ the modified counting rule Q ~ {p 7 Qq 1 , r^ 1 } for the 1 Sq channel as well. 
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and it is related to the S'-matrix via 

S -l = e 2 * - 1 = ™E = i (E^L) As , (43) 

pcot<5 — ip \ 2n J 

where So is the phase shift for the 1 S'o channel. Meanwhile effective range expansion 
(ERE) reads 

X x 

pcotS = + + ( 44 ) 

a z 

where ao (= —23.71 fm) is the scattering length and r (= 2.73 fm) the effective range in 
the 1 S'o channel. Inserting Eq. (44) into Eq. (43), one obtains a s = —1, and 

2/27 m N r 1 

V* = \ — , D s {p) = — r —. j . 45 

m N V r 2 — + ip - ^r p 2 

Now we fix the LEC's in the coupled channel. Diagrams for the 3 S 1 — 3 Di channel are 
depicted in Figs. 11 (a, b). Each of the amplitudes reads 

4tt 1 

V2 cfV i .... 

6 VW m Npdm N -_r \i ^iP ip 
The relation between the S'-matrix and the amplitudes in the coupled channel is given by 

r< .pm N / A[ss] A [DS ]\_f e 2l5 "°cos2ei te^+^hi^h \ 

b ~ 1+1 2vr U[5D] A [DD] ) ~ U^+^sh^ e^ca^ J' ^ 

where we have employed a convention for the phase shifts defined in Ref. [47]. 

Since it is known that Z\ is numerically small, we put cos(2ei) ~ 1. Around the 
deuteron pole, ERE reads [23], 

p cot5 = -7 + \pd{p 2 + 7 2 ) + • • • , (49) 

with 7 _1 = 4.319 fm and p d = 1.764 fm. Following the same steps that lead to Eq. (45), 
one obtains o t = —1, and 

2 /2tt m N p d 1 Z d 

Vt = \—, D t {p) = — — — = — - — H , 50 

where is the wave function normalization factor of the deuteron around the pole of the 
deuteron binding energy B, and the ellipsis in Eq. (50) denotes corrections that are finite 
or vanish at E — —B. Thus one obtains [19] 

Z t =T^— (51) 
1 ~ Pdl 
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Next we renormalize the coupled channel amplitude A[ds] & t the deuteron pole p = 27 
through the relation, 



_ 2 ( Vsd \ tan(2ei' 



i-VsdJ sm(5 - 5 2 ) 



2A 



[SD] 



p=il A [SS ] - A[ DD] 



-■^T 2 , (52) 



71 



where we have neglected the D-wave amplitude A\dd\- Thus, neglecting the rjl d correction 
in the relation above, one obtains 



Cf*) ~ 6,/—^. (53) 
V m N 7 
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